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The generation of electricity using direct electrostatic and magnetohydrodynamic (MHD) con- 
version of the plasma particle energy of small to mid-size chemically assisted microwave or glow 
discharge plasmas (CA-plasma) power sources in the range of a few hundred Watts to several 10's 
of kW for microdistributed commercial applications (e.g. household, automotive, light industry, 
and space based power) is studied for the first time. In the determination of the effect of plasma 
parameters on conversion efficiency, careful attention was paid to the unique plasma conditions of 
low pressure, low ionization fraction, and nonthermal ion energies that are much greater than that 
of the thermal ions of traditional MHD but much lower than those of a fully ionized plasma typically 
generated for fusion experiments. The density of plasma ions and neutrals and their cross sections 
for processes such as charge exchange were also considered. The most important parameters were 
found to be charged particle density and energy, as well as the large inventory of neutral gas atoms 
and molecules. Momentum and charge exchange of plasma ions with the large background fraction 
of neutrals represents a limitation to conversion efficiency. Two conversion technologies were ex- 
amined in some detail. We considered the possibility of converting a CA-plasma using adaptations 
of a member of the broad category of electromagnetic direct converters previously developed for 
recovery and conversion of the high energy particles lost from tandem mirror and magnetically 
confined plasmas, and an MHD converter previously developed for conversion of high pressure 
combustion gases to electricity. While it was found that both conversion techniques performed 
well under ideal conditions for conversion of plasma to electricity showing conversion efficiencies 
of > 50%, the tight coupling of plasma cell and converter, size limitations, particle energy, and 
the substantial inventory of relatively low energy neutrals eliminate direct electrostatic converters 
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as practical converters under these conditions. However, MHD conversion of CA-plasmas appears 
feasible at ~50% efficiency with a simple compact design. 
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I. INTRODUCTION 

Central station generation and distribution, the mainstay of electrical power production for the 
last 100 years worldwide, is now being supplemented in an increasing number of areas by smaller 
power units closer to the end-user group. Most distributed-generation units are in the capacity 
range of 100 kW-3 MW (electric), but some could be as large as 250 MW (electric). Distributed 
generation solves some of centralized power's inherent problems of transmission and distribution line 
losses, electromagnetic pollution fears from high-tension lines, cost and difficulty of transmission- 
line maintenance, and inefficiencies in load factor design of power plants (wherein the use of a 20% 
capacity safety factor is still a common industry practice when estimating peak loading). 

The microdistributed market emerged with uninterruptable power supplies (UPS) which serve 
the premium power market including businesses where brief electrical outages can cause severe 
monetary loss: telecommunication sites, computer centers, server hotels, e-commerce centers, semi- 
conductor fabrication facilities, and others. With a suitable technology, market conditions exist to 
extend the trend away from central power with the expansion of microdistributed power into the 
broader electricity market. The broader market which includes hundreds of millions of homes and 
businesses in the US, Europe, and Japan will be drawn by significant cost savings and increasing 
unreliability of the grid with a lack of viable microdistributed alternatives. The market is already 
moving in that direction. Plug Power 1 is developing a home fuel cell unit for the U.S. market, 
and in Japan, Matsushita plans to introduce a 1 kW fuel-cell system costing $10,000 where the 
price of electricity is 20 /kW-hr for microdistributed power. 2 The current price for government 
subsidized green distributed power in Germany is even higher, 50 /kW hr. 3 The populace of the 
third world, particularly Asia, represents a further enormous market opportunity for a technology 
with low capital and O&M costs and no requirement of an electrical grid infrastructure. 

The automotive market is also currently at a crossroads with many different options being con- 
sidered for the next generation of automobile. Evidence of the changing landscape for automobiles 
can be found in the recent increase in research into the next generation of automobiles. But, the 
fact that there is no clear front-runner in the technological race to replace the internal combustion 
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(IC) engine can be attested to by the divergent approaches taken by the major automobile com- 
panies. Programs include various approaches to hybrid vehicles, alternative fueled vehicles such as 
dual- fired engines that can run on gasoline or compressed natural gas, and a natural gas- fired en- 
gine. Serious efforts are also being put into a number of alternative fuels such as ethanol, methanol, 
propane, and reformulated gasoline. To date, the most favored approach is an electric vehicle based 
on fuel cell technology or advanced battery technology such as sodium nickel chloride, nickel-metal 
hydride, and lithium-ion batteries. 4 Although billions of dollars are being spent to develop an al- 
ternative to the IC engine, there is no technology in sight that can match the performace of an IC 
engine system. 5 

A chemically generated or assisted plasma (CA-plasma) as a novel power source using hydrogen 
as the fuel has been reported previously. 6 " 11 Since the power is in the form of a plasma, high- 
efficiency, low-cost direct energy conversion may be possible, thus, avoiding a heat engine such as 
a turbine 12,13 or a reformer- fuel-cell system. Significantly lower capital costs and lower commercial 
operating costs than that of any known competing energy source are anticipated. 

High temperature plasmas possess a substantial inventory of energy stored in the thermal and/or 
kinetic components of plasma ions, electrons, and in some cases neutral gas particles in some weakly 
ionized plasmas. There is obvious incentive then in devising methods and technologies to efficiently 
extract this energy and convert it to a more useful form. Most often, conversion to electrical energy 
is desired as this form is readily stored and transmitted, and is efficiently converted to mechanical 
work at the delivery site. 

A number of conversion schemes have been studied in the four plus decades of controlled 
thermonuclear fusion research. At high temperature (as that produced in the blanket material of 
high power D-T fusion reactor) a thermal steam cycle 14,15 is usually considered the most practical 
energy extraction means as the bulk (80%) of the energy release is in the form of chargless neutrons. 
Thermal steam cycles are robust, reliable, proven technologies, and are well established as the work 
horse of modern electrical power delivery. Yet, the conversion efficiency is limited and high coolant 
temperatures are required. Furthermore, costs are prohibitive for the use of steam cycles in small, 
distributed power sources. 
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Direct conversion of plasma charged particle kinetic to electric energy 16 may represent an at- 
tractive alternative to the steam cycle for at least several plasma systems of great interest including; 
(a) the D-T fusion reactor (as a "topping" unit to extract the 20% of fusion energy in high en- 
ergy charged particles), (b) advanced, a-neutronic fueled fusion reactors, and (c) C A- plasma cells. 
In fusion reactors, the fully ionized, high temperature (up to 10-15 keV) plasma energy may be 
readily extracted by direct, electrostatic means, thereby converting charged particle kinetic energy 
to electrostatic potential energy via decelerating electrodes. 16 Whereas for CA-plasma cell devices, 
possessing only weakly ionized and relatively cold plasmas, conversion methods more compatible 
with a fluid environment like MHD converters may be required to extract stored energy. 

A number and variety of direct energy conversion techniques have been studied over the years. 16 
Many of these may be loosely grouped into the following broad catagories. (1) Electrostatic Direct 
Converters : Electrostatic direct devices convert directed ion kinetic energy to electrical potential 
energy via an electrode (or set of electrodes) electrically biased to decelerate ions extant from 
the plasma source. The most well studied of such converter devices are the Venetian blind 17 * iS and 
periodic focused 19 converters. These devices appear to hold great promise as very efficient (80-90%) 
direct converters for large scale (on the order of 1000 MW) generating stations. In these devices, 
plasma particles are electrostatically separated before deceleration and collection at the electrodes. 
Separation incurs space charge limitations which are particularly troublesome for all but very high 
energy particles. Reasonable currents are achievable only at very high energy (several to 100s of 
keV). Particles of these energy levels are not present in appreciable numbers in CA-plasma cells. 
Furthermore, to mitigate the effects of high heat loading, 20 ' 21 such devices require plasma expansion 
and become enormous in linear scale (10's-100's of meters). 

(2) Electromagnetic Direct (Crossed Field or Drift) Converters : The guiding center drift of charged 
particles in magnetic and crossed electric fields may be exploited to separate and collect charge 
without the necessity to do so electrostatically. Space charge complications are to a large degree 
eliminated. Dimensions can often be reduced (for low power converters) by many orders (perhaps 
to the ion gyro-scale). Natural mating of the converter magnetic field to a guide field is a further 
advantage. As the devices extract particle energy perpendicular to the magnetic field, expansion 
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may not be necessary and is often undesirable. The performance characteristics of an idealized 
E x B converter which relies on the inertial difference between ions and electrons, is analyzed in 
sec. III. Timofeev^ devised a high effidency ^ ^ ^ ^ g ^ ^ 

VB drift collection. This particular device is again designed for high-power fusion energy conver- 
sion and is quite large in dimension, and requires expansion and end plug fields to prevent plasma 
leakage. In addition, collisions among energetic ions and neutral particles in a low power p.asma 
cell will likely interrupt the drift trajectory required for efficient energy conversion. As an example, 
Ar* ions in a 1 T field have a gyro-frequency of ^ ~ 2.4 x 10* B -i ^ a frequency 
with neutral Ar atoms of u in ~ 6 x 10* .-» at 40 eV, making the ion magnetization parameter 
ft = «*/* ~ 0.04. (For H+ ions under the same conditions the magnetization parameter is 0.27). 
Ions then are readily interrupted in their drift trajectory and will not reach the desired coUection 
electrode in the VB x B direction. 

(3) MHDCgnverter: I n MHD a highly conducting plasma, flowing at velocity « in a direction 
across a magnetic field B gives rise to an electric field E = -* x B. This electric field may be 
intercepted at the boundary of a plasma device by electrodes and exploited to drive electric current 
through an external load. Mechanical flow energy of the conducting fluid is then converted to 
electrical energy. In the presence of a load to complete the circuit, the density of electric current, 
it >s given by the plasma Ohm's law 

] = a(B + uxB) 

where a is the plasma electrical conductivity. The term is « x B is referred to as the MHD electric 
field or MHD term in Ohm's law. 

The performance of an MHD power conversion system is impacted strongly by the value of 
attained in the plasma region of the MHD converter. As such, collisions among charge carriers 
or between charge carriers and neutral gas atoms play a crucial role. Collisions, however, are 
not as disruptive in MHD converters as they are in direct conversion. MHD converters operate 
on fluid plasmas where collisions are frequent and the trajectories of individual plasma particles 
are relatively unimportant. The conductivity is also affected strongly by the strength of applied 
magnetic field in plasma. A more detailed discussion of this subject is presented in sec. IV. MHD 
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Power conversion devices acting on alkali metal seeded, high-pressure (at or above atmospheric) 
combustion gases have been extensively studied 24 " 29 including MHD concepts to deliver AC power 
directly from the converter. 30 ' 31 Prior studies also include plasmas generated in shock tubes or arc 
jets for high power electric generation. Little effort has been made in studying MHD conversion in 
hotter, more tenuous plasmas for low power applications. The work described herein is intended 
to initiate such a dialogue. 

II. PROPERTIES OF CA-PLASMA CELLS & CONVERTER CONCERNS 
A. Plasma Parameters 

Plasma cells operating at about 1 Torr total pressure in the glow or microwave discharge regime 
producing output power significantly greater than the input power due to catalytic release of energy 
from supplied hydrogen have been reported previously. 6 " 11 Results from studies on ca-plasmas of 
hydrogen with strontium, argon or helium with 3-10% hydrogen and strontium with these latter 
mixtures include: 

1. the observation that glow discharge plasmas of the catalyst-hydrogen mixtures of strontium- 
hydrogen, helium-hydrogen, argon-hydrogen, strontium- helium- hydrogen, and strontium-argon- 
hydrogen showed significant Balmer-a line broadening corresponding to an average hydrogen 
atom temperature of 25-45 eV; whereas, plasmas of the noncatalyst- hydrogen mixtures of 
pure hydrogen, krypton-hydrogen, xenon-hydrogen, and magnesium- hydrogen showed no ex- 
cessive broadening corresponding to an average hydrogen atom temperature of ~ 3 eV, 10,11 

2. the observation that microwave helium-hydrogen and argon-hydrogen plasmas having cata- 
lyst Ar + or He 2+ showed extraordinary Balmer-a line broadening due to hydrogen catalysis 
corresponding to an average hydrogen atom temperature of 110-130 eV and 180-210 eV, 
respectively; whereas, plasmas of pure hydrogen, neon-hydrogen, krypton-hydrogen, and 
xenon-hydrogen showed no excessive broadening corresponding to an average hydrogen atom 
temperature of ^ 3 eV, 7,10 
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3. the observation that microwave helium-hydrogen and argon-hydrogen plasmas showed av- 
erage electron temperatures that were high, 28,000 K and 11,600 K, respectively; whereas, 
the corresponding temperatures of helium and argon alone were only 6800 K and 4800 K, 
respectively, 7,10 

4. the observation that the optically measured output power of gas cells for power supplied to 
the glow discharge increased by over two orders of magnitude depending on the presence 
of less than 1% partial pressure of certain catalysts in hydrogen gas or argon-hydrogen gas 
matures, and an excess thermal balance of 42 W was measured for the 97% argon and 3% 
hydrogen mixture versus argon plasma alone, 6 

5. the observation that the power output exceeded the power supplied to a hydrogen glow 
discharge plasmas by 35-184 W depending on the presence of catalysts helium or argon and 
less than 1% partial pressure of strontium metal in noble gas-hydrogen mixtures; whereas, 
the chemically similar noncatalyst krypton had no effect on the power balance," 

6. the Calvet calorimetry measurement of an energy balance of over -151,000 kJ/mole-H 2 with 
the addition of 3% hydrogen to a plasma of argon having the catalyst Ar+ compared to the 
enthalpy of combustion of hydrogen of -241.8 kJ/mole-H 2 ; whereas, under identical conditions 
no change in the Calvet voltage was observed when hydrogen was added to a plasma of 
noncatalyst krypton, 8 

7. the observation that upon the addition of 10% hydrogen to a helium microwave plasma 
maintained with a constant microwave input power of 40 W, the thermal output power was 
measured to be at least 400 W corresponding to a reactor temperature rise from room temper- 
ature to 1200 °C within 150 seconds, a power density of 40 MW/m 3 , and an energy balance 
of at least -5 x 10 5 kJ/moIe-H 2 compared to the enthalpy of combustion of hydrogen of 
-241.8 kJ/mole-H 2> 9 

Plasma parameters of a CA-plasma of an argon-hydrogen mixture (97/3%) which are very 
conservative in terms of their impact on the potential for conversion of plasma power to electricity 
are given in table 1 along with fill gas parameters at ~ 1 Tbrr. 
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B. Device Scale 

Practical design parameters were established for the magnetic field strength B of 1 T and 
a converter physical scale L of 1 m to determine the performance of the converter. Fields on 
this order are readily produced with Weiss electromagnets 32 with iron or rare-earth cores and 
without active cooling. The choice of converter length is considered a reasonable upper limit in an 
initial analysis of micro-distributed power devices without a detailed cost analysis. At times, the 
parameters B and L will be treated as independent variables for the purpose of optimization or 
parameterization; although, it is always recognized that practical considerations limit these to be 
of the order prescribed above. 

C. Energy Content and Power Flow 

To gauge the power scales involved, initial energy content and power flow estimates were made. 
At T e ~ 10 eV, T{ ~ 40 eV t and n e? j = 10 12 cm" 3, for a 1 liter plasma ceD, the stored thermal 
energy is 

U T = n Cii V{kT € + kTi) ~ 8mJ (2) 

independent of ion species, where V is the plasma volume. Power flow can be estimated presuming 
the ability to extract this stored energy at an acoustic rate. At 40 eV, H+ ions have an acoustic 
speed of ~ 6 x 10 4 m/s, while that for Ar-f ions is ~ 10 4 m/s. For a 10 cm drift path through the 
cell, this makes the drift time for H + and Ar + ions t ~ 1.67/is and ~ IOjis, respectively. Assuming 
a particle replacement to maintain steady conditions of nvA y where nv is the particle outflux and 
A is the flow channel cross section, thermal energy can then be extracted at a rate Ur/t ~ 4.78 kW 
for H + and 0.8 kW for Ar + . This analysis is identical to setting the power output equal to the 
kinetic energy flow rate, mv 2 (nvA)/2. A similar argument can be made by considering the rate at 
which work is done on the fluid as it is expelled from the cell, P = dW/dt ^ Fv for constant F. For 
F ~ pA y then P — pvA = nkTvA for an ideal gas. All of the aforementioned arguments, of course, 
require that steady conditions are maintained while extracting energy at the estimated rate. To 
do this, plasma particles must be replaced at the rate nvA and heated to steady temperature T at 
the rate kT/t per particle. 
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D. Converter Concerns 

It is considered convention in the direct conversion of plasma to electrical energy that a conver- 
sion system must perform the following tasks. (1) A well defined plasma stream is extracted from 
the plasma confinement or reactor chamber. (2) Neutral particles must be trapped or otherwise 
diverted from the plasma fiow to ensure high quality plasma and reduce the deleterious effects of 
neutral interactions including elastic scattering and charge exchange recombination. (3) The extant 
plasma stream must be expanded to reduce the heat loading on converter surfaces such as elec- 
trodes and grids, and convert plasma thermal energy to flow energy, thereby enhancing converter 
performance. (4) Charge separation is performed usually by electrostatic means. Ions and electrons 
are separated very early in the converter region. Before being individually collected, substantial 
space charge is developed which can severely limit performance especially at low energy (below a 
few keV). (5) Charged particles are decelerated at high voltage electrodes and collected . For high 
efficiency, many electrodes may be required, each set at a different bias voltage to intercept ions 
with kinetic energy nearly equal to the bias potential. (6) To meet the needs of the application, 
direct converter power (usually high voltage DC) must be conditioned to the requirements at the 
delivery site. 

These are all challenging engineering issues for direct conversion. Neutral trapping incurs 
the need for diverting the plasma flow and a differential pumping system to remove the neutral 
inventory. Extraction may require a separate extraction chamber and additional magnetic field 
coils for guide fields. Expansion increases the physical dimensions and cost of the converter and 
supporting systems. Charge separation introduces the inevitable and sometimes fatal complication 
of space charge limitations. 

Fortunately, however, many of these requirements are dictated by the high power and high 
energies per particle associated with the nuclear fusion origins of direct conversion. At lower particle 
energies, many of these constraints can be relaxed. For example, wall loading is not considered 
a materials concern at the relatively low power typical of CA-plasma cell devices. Therefore, 
expansion may not be necessary unless there is a compelling conversion advantage (i.e. greatly 
increased conversion efficiency). Collecting plasma thermal or flow energy directly (as in the E x B 
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or MHD conversion techniques) rather than first converting this energy to directed individual 
particle energy, eliminates the need for charge separation. This represents an enormous advantage 
and, in many cases, an enabling condition for C A- plasma power converters. Space charge would 
otherwise pose an insurmountable obstacle. As an example, consider an infinitely wide (so that 
transverse space charge effects are neglected) beam of H + ions at 10 eV and 1.6 kA/m 2 . This beam 
dictates a longitudinal space charge limitation on the maximum collector length of ~ 0.5 mm, an 
obviously impractical constraint. The complications associated with plasma extraction and neutral 
trapping may also be eliminated by considering conversion strategies that allow the immersion of 
electrodes in-situ or allow a collector region to be closely coupled to the plasma cell to provide 
a natural flow path from cell to collector as in the case of E x B and MHD converters. Flow or 
collection interruption by neutral particles must still be considered, however. 

E. Recombination 

In low temperature plasmas, the recombination of free charge is often an important consider- 
ation in determining the concentration and distribution of charge states. This is especially true 
of low temperature plasmas that possess high neutral concentration and low ionization fraction. 
Three principle reactions dominate in the parameter range of interest, radiative and dielectronic 
recombination, and resonant charge exchange (CX). Estimating the rate of dielectric and radiative 
recombination is important for any plasma- to-electric conversion scheme since these processes re- 
move free charge from the inventory intended for collection. Though the CX reaction does not alter 
the net concentration of ions in the plasma, it has the deleterious effect for conversion technologies 
relying on flow like MHD by removing energetic ions from the flow stream and replacing them with 
energetic neutral particles, leaving behind cold ions. 

CX may occur among particles of the same species or among different species so long as the 
ionization energetics permits. In CA-plasma cells comprising Ar-H mixtures at low pressure in 
which the H minority concentration is not insignificant (> 3%), H + is expected to be the majority 
charge carrier as it is energetically favorable for ionized Ar to ionize H through CX. 
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The relevant recombination reactions then are those involving H+. Radiative recombination in 
H+ occurs with a rate coefficient of ~ 10" cm*/*. In plasma of electron density of 10" cm"', this 
yields a recombination frequency of 0.1 .-I or a 10 s recombination time. The recombination mean 
free path at 40 eV is then in excess of 600 km, a completely negligible process. CX, on the other 
hand, occurs at a much higher rate. The CX cross section for H+ on H» is ~ 10"" cm 2 . At 5% 
neutral H concentration, the CX mean path is then on the order of 10 cm. This is a reasonable 
dimension for the length of a small scale converter for a CA-plasma power cell. 

III. E x B DIRECT CONVERTER 

To illustrate an example of direct conversion for small to mid-scale power applications, the 
kinematics expressions and conversion efficiency are derived for the simple E x B converter with 
both ion and electron collectors. The zero order behavior of an ideal converter is described to 
retain an analytically tractable formalism. In the absence of significant expansion, the impact of 
collisions cannot be underestimated. As briefly mentioned in the introduction section, collisions may 
significantly reduce the efficiency of such devices by interrupting ion trajectories to the collector. 

A schematic of a converter based on E x B collection is shown in Fig. 1. Here a rectangular 
arrangement of electrodes is chosen for simplicity with plasma particles incident from the left 
and drifting along guide field, B. When both ions and electrons enter the collector region and 
experience the applied crossed fields, g and B, they will immediately assume a guiding center drift 
in the direction perpendicular to both E and B and with speed v B = ExB. Though this speed 
is identical for ions and electrons, ions having greatly reduced transitional speed parallel toB (for 
Ti = T e ), will be turned and deflected to the upper electrode before electrons. For the same transit 
time, high-speed electrons will then be intercepted by the end electrode. The addition of electron 
collection increases direct power conversion. 

Direct electromagnetic conversion like E x B offers distinct advantage over electrostatic con- 
version for a number of reasons. Perhaps the single most important aspect is that, like all fluid drift 
conversion processes, E x g conversion acts on the entire neutral plasma. The necessity to separate 
charge is thereby removed as are the space-charge complications that arise therefrom. Coupling to 
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where M is the species mass. An expander region must conserve magnetic flux so that for cross 
sections of linear dimension * and „ 2 at the inlet ^ M ^ ^ ^ 

B>d*=B 2 4 (4) 
Assigning * = adj> where a is a dimensionlesS) inverse expansion ^ ^ ^^.^ ^ ^ 
expressed as S 2 = B y conserving ^ ^ = ^ ^ ^ 

post expansion particle speed are obtained 

Note that when a = 1 (the no expansion limit) then w Xl = v±l and tr,,, = ^. 

Conserving mass flow „„„ (where -4 is the channel cross section at any position) determines 
the particle density change across the expander section 



"2 O 2 



"i (3 - 2o 2 )V2 (6) 

Coupling this result to an adiabatic expansion requirement {pin =const.) and the ideal gas law 
(P = nkT), dictates a temperature difference across the expander region 

l^ 3 -^^) 7 ^- 2 (7) 
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where h and / are the length* of the ceU and expander regions, restively. At a = 1/2, h/l = 1/2, 
and 7 = 5/3, the relative temperature decrease is found to be T 2 /T x ~ 0.2. 

Ion energy extraction requires that ions drift to the ion collector at E x B. The rate of ion 

collection is then 



where A ~ 2d 2 is the width of the collector electrodes. Since ions bring only perpendicular energy 
W u to collection, the rate of energy collection is 



Pi = n 2 v ExB lAW,, = — n ,kT, —I A. 

3 (3-2o2)»/2 ni * Tl ^ ,A (9) 

Electrons by contrast bring W h to the electron collector since they drift parallel to B. Ambipol 
considerations require electrons to reach their collector at the same rate that ions reach the 
collector, so that 

*2tWA = mvn u & 2 (10) 
so that the electron drift speed is limited to ^ = (//A)*,,,. The collected electron power is then 

^ \2 mV ^J (3-2a*)V2 2 aUJ (") 
Combining ion and electron power, the total collected power is 

For l/dx » 1, P peaks at a „ ^ Since ^ fa ^ ^ ^ 

1, expansion is not desirable in this situation. At • - l, there is no adiabatic change in fluid 
parameters across the expander region, and 

separating the ion and electron contributions, respectively, within the square brackets. 

The power input to the converter has two contributions. Thermal flow power from the cell is 

estimated 

P <<™ = 3nkT \f%~ (14) 
14 
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and in addition, there is a contribution from the work required of an external agent (power supply) 
to maintain the electric field in the converter in the presence of particle drifts 



By the mass ratio M 3> m, the ion contribution is the only important contribution to this compo- 
nent of input power. 

The converter efficiency is defined r) = P/Pm where Pj n - P/iaw + Pe- Defining a new param- 
eter as the dimensionless drift- to- thermal-speed-ratio, a = v BxB /ve^, the efficiency expression is 
parameterized as 



With l/d ~ 5 the conversion efficiency peaks at r) ~ 70% near a *>* 1. This is demonstrated in 
Fig. 2 as a parameterization of rj vs. a with l/d~ 5. 

Some additional conditions should be considered. First, to avoid transverse ion loss, it is 
necessary to ensure rj,. < A/2 — the ion gyro-scale fits within the channel dimensions. This 
requires B > V2kT/qd ~ 45 G for 10 eV H + ions at d = 10 cm, a trivial requirement. More limiting 
is ensuring a drift time much larger than the ion gyro- time to allow fully developed ion drift flow 
to intercept the upper electrode (&VB/E) > u/" 1 . This places an upper limit on E < &qB 2 /M of 
£,7600 V/m when B = 200 G. Since l/d > 1 is required for V|j Je > v SxS , this forces f ~ 1/2-1 m. 
The condition for equal ion and electron contributions to output power is 



again for 10 eV and 200 G. This is superceeded by the gyro-time requirement. It is more reasonable 
then to fix the electric field to a lower value near a ~ 1. For example, at E = 1000 V/m (e.g. V =. 
200 V for A = 20 cm) and B = 200 G at 10 eV, one finds a - 1.6 and q - 54% (from Eq. (16)). 
The power output under these conditions (and with n t ~ 10 12 cm" 3 ), is P ~ 4.7 kW. 

A reasonable quantity of electric power may be extracted in such a converter design provided, 
of course, that collisional effects are not important. The effect of collisions has not been considered 
here. One may estimate, however, as done in the introduction (sec. I) that the magnetization 




(15) 




(16) 




(17) 
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parameter is quite low for ions, 0 f 1. This is found even at very high field approaching 1 T. 
The majority of ions in the stream are then not expected to complete an uninterrupted trajectory 
to the ion collector. For this reason alone, the success of direct collection in a compact, low power 
converter would be suspect. 



IV. MHD CONVERTER 



The MHD converter exploits the Lorentz action on a flowing magnetofluid across a magnetic 
field to generate an electrical potential difference. A schematic illustrating the basic components 



of an MHD system interfaced with a CA-plasma cell device is shown in Fig. 3. Magnetofluid flow 
is extant from the CA-plasma cell (labeled Plasma Tube-Reactor) at flow velocity u. As the flow 
enters the MDH converter/expansion region, it experiences flow-crossed magnetic field B. In the 
absence of an external load, an open-circuit electric field E Q - -ux B is generated. This is a 
direct expression of the plasma Ohm's Law when the flow of electric current is prevented. When 
the system supports the flow of electric current, the relationship between electric current density / 
and E (Ohm's law) is given in Eq. (1) where a is the electrical conductivity of the magnetofluid. 
The circuit is completed through the external load which reduces the electrode voltage so that 
E = kE q = -ku x B where k < 1, so that the magnitude of the current density becomes 



The continuous appearance of an MHD voltage, V = Ed, (where d is the electrode separation) 
and electric current flow is predicated upon continuous fluid flow through the channel defined by 
the separation of the converter electrodes. The flow may be maintained via a pressure drop, VP, 
across the channel so that the plasma component of the fluid is in dynamic equilibrium with the 
applied field 



or Ap = jBL for a linear pressure drop (i.e. constant B,j) across a channel of length L. The 



j = <r(l - k)uB 



(18) 



Vp = fxB 



(19) 



converter length required to support the fluid at fixed Ap can then be written 




(20) 
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indicating a reduction in scale for concomitant increases in u, B, or a. When the field magnitude 
and flow speed are fixed by power and materials limitations, there is a premium on a high degree 
of conductivity for the fluid. The electrical conductivity in plasma is determined by a number of 
factors including species, charge, average thermal speed, collision cross section, 33 and B. A detailed 
analysis on conduction in partially ionized gases will follow (sec. A). 

In order to maintain a pressure drop, and hence flow, evacuation of the fluid extant from 
the converter is required in conventional MHD, but may be unnecessary in the CA-plasma case. 
Three possible scenarios are presented. (1) The plasma cell and converter may be directly coupled 
and open to atmosphere in a once-through, open system. The pressure drop is maintained by a 
vacuum pump or by operating at greater than atmospheric pressure. (2) The cell and converter 
may be arranged in a closed configuration which utilizes a recirculating pump to accumulate the 
converter effluent and divert it back to an injection reservoir in the plasma cell. Neither of the pump 
scenarios presented in (1) or (2) are beneficial in an energy conversion system since the pumping 
power required to maintain Ap and hence u would necessarily be greater than that converted to 
electrical power by the flow. This follows directly from the requirement that the energy extracted 
as electricity is a direct consequence of a high enthalpy fluid expending flow energy in crossing 
B. (3) Hot plasma generated in the CA-plasma cell and expanding outward therefrom into the 
converter region, may introduce an adverse pressure gradient which may be filled by backflow of 
neutral gas from the converter region returning to the cell. A natural convection like pattern may 
be established providing both continuous flow and refueling simultaneously. The CA-plasma cell 
and converter may then be coupled in a simply closed configuration without need for pumping. 

Power flow through an external load at MHD supported / and E can be computed 

P = j-E = ck{1-k)u 2 B 2 (21) 

This is optimized at k = 1/2, which simply represents the impedance matching condition where 
half the open-circuit voltage drop appears across the load. 
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A. Electrical Conductivity in an Applied Magnetic Field 

As the electrical performance of the MHD converter is a strong function of the electrical 
ductivity, the accuracy of quantitative determination of a is of paramount importance. Indeed, 
when the high concentration of neutral particles is in flow and thermal equilibrium with plasma 
ions and electrons, it is noted" that the MHD efficiency of conversion in historical applications is 
limited to the ionization fraction. This is a rather debilitating limitation as the ion fraction may 
be quite low, perhaps at only a few percent or less. However, in the low pressure CA-plasma case, 
no such equilibrium exists. The greatly reduced collisionality afforded by low density somewhat 
decouples the plasma species from neutral particles. Input power then is not required to heat the 
large inventory of neutrals, nor is it required to drive flow in this component so that the input 
power requirements may be much reduced and the electrical efficiency much greater. 

The strong applied magnetic field, on the other hand, does have a dramatic influence on con- 
duction.^ This can be directly ascertained from the electron momentum equation. Under the 
limiting, yet illustrative, conditions of constant and uniform plasma and flow with B = Bz and 
VT e = 0, the electron equation of motion yeilds the familiar expression for the transverse electron 
speed 

1 + ftfB ^l+figeB n e " TTflf^ ~ 1 + Og (22) 
where p. = e/(m e u e ) is the electron mobility, D e = W e /(m e „ e ) is the electron mass diffusivity, and 
ft e = o/c. = eBI(m e u e ) = n e B is the electron magnetization parameter. For ft e » 1, the electron 
fluid is magnetized and strongly influenced by B. When fl e « 1, the applied field has much less 
influence than particle collisions on electron transport. The first two terms in Eq. (22) are the 
familiar electric and diamagnetic drift terms perpendicular to B. The last two terms represent 
electrostatic mobility and diffusive transport, yet the magnitude of the transport coefficients is 
reduced by the factor 1 + 0*. so that 

^-ITnr D " = TTm (23) 
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which may represent a significant reduction for large B. Alternatively, the effective collisionality is 
increased 

u €L = u t {\ + (l 2 e ) , q eL = r? e (l + flj) (24) 

Because of the mass difference, ions and electrons are influenced by collisions and B to differing 
degrees. Table 2 shows ion and electron collision frequencies with all species present (e, H+ , Ar° 
neutrals) for a low power plasma cell possessing conditions of table 1 with 40 eV ions, 10 12 cm -3 
plasma density, and a Coulomb logarithm of 20. Charge neutral collisions are among ions or 
electrons with neutral Ar atoms at 1 Tbrr total pressure. Coulomb collisions are self (i-i or e-e) or 
cross (t-c or e-t) involving electrons and H + ions as the dominant ionized species. Conduction for 
both ions and electrons is limited by collisions with neutral particles due principally to the large 
inventory of neutrals at 1 Tbrr. These mechanisms (i.e. i/ cn , v in ) will then be considered the only 
important collisional effects. 

In weak magnetic fields (fi t \e 1)» both ions and electrons are relatively unaffected by the 
presence of B. Electrons, then, due to their higher mobility, dominate electrical conduction. When 
the magnetic field strength increases such that n tt€ » 1 (the strong field limit), both ions and 
electrons are magnetized and electrical conduction perpendicular to B is dominated by ion flow. 
For intermediate fields, as for our test case near B ~ 1 T, both ions and electrons conduct electrical 
current. To see this, a conductivity ratio can be estimated 

f*.2#±<S> (25) 

At 1 T, fi e ~ 29.3 and ft, ~ 0.27 (H + ions) while o^Oi = ~ 85.7. The conduction ratio 

perpendicular to B then becomes o t Joi L ~ 0.1 so that electrons carry only about 10% of the 
electrical current in the converter. 

B. Sample MHD Converter Performance 

MHD converter performance can be illustrated by examining a test case with B^1T,k = 1/2, 
Ap = 1 Tbrr, and u = 1.36 x 10 4 m/s (ion acoustic speed at neutral Ar inertia). For v CL = 0.1ft x , 
and n e = 10 12 cm -3 , the plasma conductivity perpendicular to B is estimated at 0.048 mho/m or 
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r?x ~ 21 n-m. Then employing Eq. (20), a converter length of only ~ 40 cm is found. This is a 
modest requirement and suggests that such large fields may not be necessary. The MHD electric 
field generated in this case is E = kuB ~ 6.8 kV/m, providing a voltage drop of 680 V across 
a 10 cm converter gap. The electric current density can then be estimated from Eq. (18) to be 
~ 326 A/m 2 so that the MHD output power becomes ~ 8.8 kW. 

Sample converter performance is shown in Figs. 4 and 5 illustrating MHD current density and 
power as functions of applied field, B, from 0-40 T at constant L and u. Though the upper limit 
on B is clearly impractical, the range encompasses all the relevant MHD physics. Under these 
constraints the MHD voltage and EMF are linear functions of B. The electric current density 
(Fig. 4), however, shows much more interesting behavior. There are two peaks in the curve, one 
at B ~ 1/ Me and another at B ~ before asymptotically decreasing to zero as B -* oo. 
This behavior is explained by considering the perpendicular conductivity or mobility of charges in 
strong B. At low field (B < l//i e ), electrons easily conduct under the influence off;, and j increases 
linearly with B since the MHD E increases with B. As B approaches l/fx e = 0.034 T, electrons 
become magnetized and are greatly impeded in their flow perpendicular to B, so that j decreases 
rapidly. This situation is sometimes referred to as the ton slip condition 27 since ions continue to slip 
through the applied field whereas electrons are trapped. In the region of B parameter space between 
\/Hc and l/fii, there is a competition between conductivity reduction and increasing EMF with B. 
At B = l/ w ~ 3.7 T, ions now become magnetized, and the current once again peaks. Beyond this 
field strength, the current is a continuously decreasing function of B. The MHD power (Fig. 5) is a 
continuously increasing function of B in spite of the variation in j with B, since E is continuously 
increasing. The power function reaches an asymptotic value (P«, = *(1 - K)<PLu 2 en/fn ~ 110 kW 
for this case) at high field since the E increase is linearly with B and ; decreases like B~ l at large 
B. 

Though the output power can approach appreciable levels, the quantity of total electric current 
remains low, / < 25 A (and ; < 600 A/m 2 ) so that induced fields remain negligible in comparison 
with the applied field. The magnetic Reynolds number Rm = ft 0 c ± ud ~ 10" 5 determines the scale 
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of flow interactions with the applied field. Since Rm <§C 1, the complications usually associated with 
flow-field distortion, hydromagnetic waves, and instabilities can be avoided. 

C. Channel Hydrodynamics 

The hydrodynamics of one dimensional, steady channel flow (u = ux ) in crossed field (B = Bz) 
is examined by considering the conservation equations of hydrodynamics 

energy : pu^ + C p tJ = y E (26) 

du - 
momentum : pu— + Vp = /xB 
ax 

continuity : m = puA = const. 

for a fluid of mass density p in a channel of cross section A. If constant flow (u = const.) is 
considered, the hydrodynamics equations simplify to 

energy : pu ± (C p T) = f-B or g = tlL (27) 
momentum : Vp = fx B 
continuity : pA = const. 

where h is the specific enthalpy. 

In considering a constant applied field and disregarding flow distortion of the applied field as 
indicated by the tiny order of the magnetic Reynolds number, the pressure profile must be linear. 
By integrating the momentum 

p(x) = 2j y £L[l-x/2I] (28) 

on 0 < x < L. Given the pressure profile above and ideal gas behavior, p = nW, the energy 
equation can be integrated to find the temperature profile 

T(x)/T 0 = [1 - x/2Lf^*^ (29) 

where T 0 is the channel inlet temperature. Since E = kuB, the exponent reduces to nk/mC p ~ 0.33 
using the properties of Ar as the bulk species. At the channel exit we find T(L)/T Q ~ 1.26. The 
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temperature increase along the channel is attributed to Joule heating of the plasma by the MHD 
current and field. 

Since the temperature and pressure profiles are determined, the density profile can be found 

"(*) = ^[l-x/2L] l+ ^ (30) 
= Ml-x/2Z,] 1 - 33 

The extant flow (at x = L) then has density reduction n(L)/n 0 ~ 0.4. 

By mass conservation, the channel cross section must widen to support constant flow while the 
density decreases, A ~ 1/p. Then it is readily determined 

«*-j&t*-*i*r***# (3.) 

= A 0 [l - x/2L)~ 1M 

As the flow exits the channel, the gap must widen to A(L)fA. ~ 2.5 or d(L)/d 0 ~ 1.6 to accom- 
modate constant flow. 

D. Generator Efficiency and the Hall Effect 

The preceding analysis considers only the idealized behavior of an MHD converter, that is the 
electrodynamic and hydrodynamic behavior in the absence of heat and particle losses, and Hall 
currents. Collisions, on the other hand, are fully accounted for through the explicit determination 
of collision frequency and its implementation in the Ohm's Law (Eq. (1)). In this context, the MHD 
efficiency may be quantified by considering the following. The output power density is determined 
by the rate at which specific enthalpy in the flow is converted to electrical energy 

dh _ - . n 
fm fa = 3 E = hE (32) 

The rate at which energy is expended is attributed to work done by the fluid in expanding through 
the applied magnetic field 

h % = n h B (33) 
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The ratio of these two expressions is the MHD efficiency 

■\ E 

Vmhd = ^ = * (34) 

This quantity is a constant (k — 1 /2) in the heretofore provided formalism since no physical effects 
other than impedance matching are considered. 

At high applied field strength, however, the Hail effect may play an important role in channel 
dynamics. 36 " 3 * The Hall effect in plasma is a consequence of electric current interaction with applied 
magnetic fields just as that experienced in metalic conductors. In plasma, though, this effect can 
have significant impact on plasma impedance and dynamics. The Hall effect is quantified via the 
generalized Ohm's Law 

E + u x B = + — (fx B - Vp e ) (35) 
a en \ / 

where the last two terms were omitted in the form previously employed (Eq. 1). The fxB term on 
the right side is the Hall term. The last term represents the electrodynamic influence of electron 
pressure gradients and is ignorable when /? e = 2^ 0 kT e /B 2 1. For C A- plasma cell conditions 
0 e ~ 10~ 6 at 1 T so that neglecting electron pressure is well justified. 

The Hall contribution, however, is most often not negligible. It can have quite a strong influence 
on plasmadynamics, electrodynamics, and energy balance in MHD plasmas. 36 " 38 Via the Hall term, 
a component of electric current and field perpendicular to B is introduced. For the cartesean MHD 
channel described earlier with B = (0,0, B) and u = (u,0,0) resulting in MHD fields E y and j VJ 
the Hall contribution appears in the -z direction and is given by 

E ha u = E X = ij x + — j y B (36) 
(7 en 

where j x is the Hall current. Lacking experimental guidance or further theoretical constraints on 
j Xi parameterization is provided by the Morozov 37,38 Hall parameter 

H = jx/enu (37) 

The MHD efficiency expression is then suitably modified to incorporate the rate at which energy 
is expended in driving Hall currents 

Vmhd = *, , ) tEx ( 38 > 
1 + uJTB 
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No credit is taken here for the potential for power conversion of the Hall current component. 
This has enhanced performance been suggested elsewhere 27 * 29 and should be considered further to 
improve the overall performance of the converter. Figure 6 displays the MHD efficiency including 
Hall losses as a function of the applied field for E = 0.1,0.3, 1.0. A smaller Hall parameter is clearly 
desirable here. As E is increased, a greater fraction of converter power is diverted to drive Hall 
currents. The effect is increased with increasing B since Ehaii ~ B for large B, 

E. Fiow 

In the absence of spontaneous plasma flow from the hot CA-plasma cell to the relatively cold 
MHD converter section, a directional plasma flow may also be formed by using a magnetic mirror. 
A magnetic mirror has a magnetic field gradient in the desired direction of ion flow where the initial 
parallel velocity of plasma particles increases as the orbital velocity decreases with conservation of 
kinetic energy and adiabatic invariant ^ = W±/B, the linear energy being drawn from that of 
orbital motion. The adiabatic invariance of flux through the orbit of an ion is a means to form a 
flow of ions along the field with the conversion of W± to 

Plasma is selectively generated in the center region of the CA-plasma power cell. A magnetic 
mirror located in the center region causes electrons and ions to be forced from a homogeneous 
distribution of velocities at the ceU center to a preferential velocity along the axis of the magnetic 
mirror. Thus, the plasma ions have a preferential velocity along the field and propagate into the 
MHD power converter. By preserving the adiabatic invariant, the parallel velocity at any position 
along the z-axis is given by 

v ll = v o-t>Lj^ (39) 
where the zero subscript represents the initial condition at the cell center. In the case that t/jj^ = 
V U ~ °* 5v ? m & J; = 0.1 at the MHD power converter, the particle velocity is 95% parallel to the 
field at the converter. 
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V. CONCLUSIONS 

Preliminary investigation has been made on plasma- to-electric conversion technologies for small 
to mid-scale C A- plasma cells. Direct electromagnetic and MHD conversion technologies have been 
considered in some detail for the unique parameter range of CA-plasma cells. The plasma conditions 
inherent in these devices (intermediate temperature and ionization fraction, high neutral inventory) 
pose unique challenges to conversion. An E x B direct converter is considered in detail showing 
quite promising ideal performance with conversion efficiency up to ~ 70%. However, in light of 
the large neutral inventory present, direct conversion including E x B, other drift concepts, as 
well as electrostatic collection should be considered as impractical as plasma-to-electric conversion 
technologies. Collisionai interruption of ion trajectories become the most serious shortcoming. 
However, fluid conversion strategies, like MHD, appear to be much better suited to cell conditions. 
These conversion technologies do not require charge separation, nor is expansion required in most, 
cases. It is discussed herein that ion acoustic flow at 40 eV can generate up to several kW of 
electrical power in a 1 T field with conversion efficiency approaching 50% including Hall losses. 
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Table 1: Typical plasma conditions in Ar-H C A- plasma cell with 5% minority H concentration. 



Electron Temperature 


T« 


10 eV 


Ion Temperature 


Ti 


30-40 eV 


Plasma Density 


n e = m 


10 12 -10 14 cm" 3 


Majority Neutral Density 


"Ar 


3 x 10 16 cm" 3 


Minority Neutral Density 




1.5 x 10 15 cm" 3 



Table 2: Ion and electron collision frequencies in a CA-plasma cell having the conditions of table 1 
with 40 eV H + ions, 10 12 cm" 3 plasma density, and Coulomb logarithm of 20. 





Neutrals 


Electrons 


Ions 




fx* s* 1 






Electrons 


6x 10 9 


2.7 x 10 6 


2.7 x 10 6 


Ions 


3.8 x 10 8 


2x 10 3 


7.8 x 10 3 
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Figures 

FIG. 1. E x B type direct converter schematic. 

FIG. 2. Ideal efficiency for an E x B converter as a function of drift to thermal ratio a with l/d ~ 5. 

FIG. 3. Schematic of an MHD converter interfaced to a CA-plasma ceil, labeled Plasma Tube- 
Reactor. A magnetic mirror coU provides a guide field for the plasma flow extant from the cell 
region. 

FIG. 4. MHD converter current density for sample converter as a function of applied magnetic field 
strength. 

FIG. 5. MHD converter power density for sample converter as a function of applied magnetic field 
strength. 

FIG. 6. MHD converter efficiency as a function of applied magnetic field strength including Hall 
losses for S = 0.1,0.3, 1.0. 
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Figure 2 
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Figure 3 
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Figure 6 
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